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Abstract 

We establish the relation between local stability of equilibria and slopes of critical curves for a specific 
class of difference equations. We then use this result to give global behavior results for nonnegative 
solutions of the system of difference equations 

b\x n 

x n +i = — ; h hi 

1 + Xn + ci y n 

n = 0, 1,..., (x ,y ) e [0, oo) x [0, oo) 

Vn+1 = — ; h h 2 

with positive parameters. In particular, we show that the system has between one and three equilibria, 
and that the number of equilibria determines global behavior as follows: if there is only one equilibrium, 
then it is globally asymptotically stable. If there are two equilibria, then one is a local attractor and the 
other one is nonhyperbolic. If there are three equilibria, then they are linearly ordered in the south-east 
ordering of the plane, and consist of a local attractor, a saddle point, and another local attractor. Fi- 
nally, we give sufficient conditions for having a unique equilibrium. 

Key words: difference equation, rational, global behavior, global attractivity, model, competitive sys- 
tem, Leslie-Gower model. 
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1 Introduction 



The study of dynamics of difference equations often requires that equilibria be calculated first, followed by 
a local stability analysis of the equilibria. This is then complemented by other considerations (existence 
of periodic points, chaotic orbits, etc.). If the analysis is applied to a class of equations dependent on one 
or more real parameters, the task is complicated by the fact that a formula is not always available for 
equilibria, and even if it is, determination of stability properties of parameter-dependent equilibria may be 
a daunting task. 

In this paper we address the study of stability character of equilibria for a parameter dependent class of 
systems of difference equations. For this we develop local geometric criteria for planar systems of difference 
equations, which under certain conditions determine local stability character of an equilibrium point. This 
idea is used to study a class of difference equations in the plane, for which we show that the number of 
equilibria determines, essentially in a unique way, the local stability character of the equilibria, as well as 
the kind of global dynamics that is possible. 

The class of difference equations we will focus on has its roots in biology applications. The Leslie-Gower 
Model in difference equations is the two-species competition model 

bi x n 

x n +i = 

1 + en x n + C12 y n 

71 = 0,1,..., (x Q ,y ) e [0, oo) x [0,oo) (LG) 

b2Vn 

Vn+1 - ~. : ; • 

1 + c 2 i y n + C 2 2 x n 

It is a modified version of the Beverton-Holt equation which is a well-known difference equation in Popu- 
lation Biology. The normalized Leslie-Gower Model 



x n+l 



Vn+1 



l + x n + ay n 

b 2 y n 



n = 0,l,..., {xo,yo) £ [0, oo) x [0, oo) (NLG) 



1 + Vn + C2 Xn 

was studied in detail by Liu and Elaydi [To], J. Cushing et. al [T], and Kulenovic and Merino [TT] . 

In this paper we consider the system of equations obtained by adding positive constants to the right- 
hand-side of the Leslie-Gower equations: 



bix n 

Xn+l = — ; + Hi 

1 + en x n + C12 y n 

Vn+1 = — ; + H 2 . 

1 + C 2 1 X n + C 22 y n 



71 = 0,1,..., (xo,yo) E [0, oo) x [0, oo) (LGI) 



The positive constants Hi and H 2 in (LGI) may account for immigration. The change of variables x n 
^jj x n , y n = ^- y n normalizes (LGI) as follows: 

bix n 

x n +i = ^— ; h hi 

1 + x n + ci y n 

hVn . , 

Vn+1 = — : + h 2 

i + y n + c 2 x n 

where in (LGIN) the tildes have been removed from the variables to simplify notation, and 

ci = — , hi = en Hi , c 2 = — , and h 2 = c 22 Hi . 
c 22 en 



n = 0,l,..., (xo,yo) G [0, oo) x [0, oo) (LGIN) 
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The system (LGIN) is an example of a competitive system, which is defined next. Let I and J be 
intervals of real numbers and let / : / — > I and g : J —* J be continuous functions. Consider the system 



•En+l — fibrillin) 

Vn+i = g{x n ,y n ) 



n = 0,1,2, ... 



(z ,yo) e/xJ 



(1) 



System ([I]) is competitive if /(x, y) is non-decreasing in s and non- increasing in y, and g(x,y) is non- 
increasing in y and non-decreasing in x. The map T(x,y) := (f(x,y),g(x,y)) associated to a competitive 
system ([1]) is said to be competitive. System (pQ) is strongly competitive if it is competitive, with strict 
coordinate-wise monotonicity of the functions f(x,y) and g(x,y). If T is differentiable on an open set 
R, a sufficient condition for T to be strongly competitive on R is that the Jacobian matrix of T at any 
(x,y) G -R has nonzero entries with sign configuration 



Competitive systems of the form ([I]) have been studied by many authors [2], [6], [13], [18] and others. The 
term competitive was introduced by Hirsch [5] (see also for systems of autonomous differential equations 
x\ = Fi(x±, . . . , x n , t)„ i = 1, . . . , n satisfying < 0, i ^ j. The main motivation for the study of these 
systems is the fact that many mathematical models in biological sciences may be classified as competitive 
(or cooperative) [16] . [T7] . [13] . Consideration of Poincare maps of these systems leads to the concept of 
competitive and cooperative in the discrete case. 

Denote with -< se the South- East partial order in the plane whose nonnegative cone is the standard 
fourth quadrant {(x,y) : x > 0, y < 0}, that is, (x±, y±) < se (x2, 2/2) if and only if x% < X2 and y± > j/2- 
Competitive maps T in the plane preserve the South-East ordering: T{x) ^ se T(y) whenever x ^ se y. The 
concept of competitive (for maps) may be defined in terms of the order preserving properties of the maps. 
Finally we note that the terms equilibrium (of a planar system of difference equations) and fixed point (of 
a map) are used interchangeably in this paper. 

This paper is organized as follows. In Section [2] a relation between local stability of equilibria and 
slopes of critical curves for a specific class of difference equations is established. In Section [3] it is shown 
that every solution to Eq.(LGIN) converges to an equilibrium. Section |4j contains the main result of this 
paper. It states that Eq.(LGIN) has between one and three equilibria, and that the number of equilibria 
determines global behavior as follows: if there is only one equilibrium, then it is globally asymptotically 
stable. If there are two equilibria, then one is a local attractor and the other one is nonhyperbolic. If there 
are three equilibria, then they are linearly ordered in the south-east ordering of the plane, and consist 
of a local attractor, a saddle point, and another local attractor. Finally, in Section [5] we give sufficient 
conditions for Eq.(LGIN) to have a unique equilibrium. 

2 A preliminary result on critical sets 

Our first result establishes a connection between local stability of the fixed point of a planar map and the 
slopes of certain curves at the fixed point. The result will be useful in the proof of the main result in 
Section [JJ 

Theorem 1 Let R be a subset 0/K 2 with nonempty interior, and let T = (/, g) : R — > R be a map of class 
C p for some p > 1. Suppose that T has a fixed point (x,y) € int R such that 




a-=fx(x,y), b:=f y {x,y), c:=g x (x,y), d:=g y (x,y) 



3 



satisfy 

< a < 1, bc>0, 0<d<l, l + (a + d)+ad-bc>0 (2) 
Let C\, C 2 be the critical sets 

Ci := {(x, y) : x = f (x, y)} and C 2 := {(x, y) : y = g (x, y)} 

Then, 

i. There exists neighborhood I C M o/x and JcM o/y suc/i t/iat i/te sets Ci n (/ x J) and C2 D (/ x J) 
are the graphs of class C p functions yi(x) and y 2 {x) for x € I. 

ii. The eigenvalues Ai, A2 of the jacobian matrix ofT at (x,y) are real, distinct, and the eigenvalue with 
larger absolute value is positive. If X\ is the larger eigenvalue, then 

- sign ( 2/1 (x) - y' 2 (x)) if b < 
- 1< A 2 < 1 and sign ( Ai - 1 ) = { (3) 



sign (y[(x) - y' 2 (x)) ifb>0 



Proof. 



i. The existence of of I and J and of smooth functions yi(x) and y 2 (x) defined in / as in the statement 
of the Theorem is guaranteed by the hypotheses and the Implicit Function Theorem. Moreover, 

yi( x ) = —r, — r~ and V2\ x ) = 1 7 — r , xei (4 

fy{x,y) l-g y {x,y) 

ii. The characteristic polynomial of the Jacobian of T, 

p(A) = A 2 -(a + d)A + (ad-6c) 

has positive discriminant, thus its roots Ai, A2 are real. Since X± + X 2 = a + d > 0, then the root Ai 
with larger absolute value is positive, and | A2 1 < Ai. Now by the hypothesis a, d £ (0,1) it follows 
that Ai + A2 = a + d<2, which implies A 2 < 1. Since p(X) has at least one positive root and since 
by hypothesis p(— 1) = l + {a + d) + ad — bc> 0, we have — 1 < A2 . To prove the second part of ([3|) , 
note that from we have 

, z_x , /_\ 1 — a c 1 — (a + d) + a d — b c 

*w-y*w = —-— d = w^d) 

(5) 

p(l) _ (A 1 -l)(l-A 2 ) 
b(l-d) b{l-d) 

The result is a direct consequence of ([5]), the inequality A2 < 1 and hypotheses (|2|). 

□ 
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3 Every solution converges to an equilibrium 

In this section we show that every solution to Eq.(LGIN) converges to an equilibrium. Let 

hx , , , , b 2 y 



f(x,y) 



+ hi, and g(x,y) 



1 + x + cxy 1 + y + c 2 x 

Then the map T(x,y) = (f(x,y),g(x,y)) associated with (LGIN) is 

hx , , b 2 y 



+ h 2 



T(x, y) 



+ hi, 



l + x + ciy ' l + y + c 2 x 
For future reference we give the jacobian matrix of T at (x,y): 

( hi (ci y + 1) 61 ci x 

Jt{x, y) = 



(x, y) G [0,oo) x [0,oo) 



(x + ci y + l) 2 
b 2 c 2 y 



(x + ci y + l) 2 
b 2 (c 2 x + 1) 



(6) 



(7) 



\ (c 2 x + y + l) 2 (c 2 x + y + l) 2 J 
By direct inspection of ([7]) we obtain the following result. 

Lemma 1 The system of difference equations (LGIN) is strongly competitive on [0, 00) x [0, 00). 
The following lemma has an easy proof which we skip. 
Lemma 2 T( [0, 00) x [0, 00) ) C [hx, hi + bi] x [h 2 , h 2 + b 2 }. 

Theorem 2 Every solution qfEq.(LGIN) converges to an equilibrium. 

Proof. It is easy matter to show that the map T is one-to-one. From ([7]) the determinant of Jx(x,y) is 

b\ b 2 {c 2 x + ciy + 1) 



det (Jacx(x, y)) 



(c 2 x + y + l) 2 (x + c\ y + l) s 



(8) 



which is clearly positive for (x, y) € [0, 00) x [0, 00). If follows that the map T satisfies hypothesis (H+) in 
|18j . By Lemma 4.3 of [18] and Theorem 4.2 of [18] we have that every solution of Eq.(LGIN) is eventually 
coordinate- wise monotone. Since by Lemma[2]every solution enters the compact set [hi, hi + 61] x [h 2 , h 2 + 
b 2 ], the conclusion follows. □ 



4 Number of equilibria and global behavior 

By solving for y and x, respectively, in the equations defining the critical sets C\ = {(x, y) G R 2 : x = 
/ (x, y)} and C 2 = {(x, y) G M 2 : y = g (x, y)} and taking derivatives we obtain the following result. 

Lemma 3 All branches of the sets 

Ci = {{x, y) G M 2 : x 2 + c\ xy + (1 - bi - hi) x - ci hi y - hi = 0} 

C 2 = {(x, y) G M 2 : y 2 + c 2 xy + (1 - b 2 - h 2 )y - c 2 h 2 x - h 2 = 0} 
are the graphs of decreasing functions of one variable. 
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Lemma 4 The map T satisfies hypotheses ([2]) of Theorem^ 



Proof. Set a := f x (x,y), b := f y (x,y), c := g x (x,y), d := g y (x,y). Implicit differentiation of the equations 
defining C\ and C 2 in ([3]) at (x, y) gives 

y'i(x) = and y' 2 (x) = -^—^ (9) 

From Lemma O y[(x) < and y' 2 (x) < 0, and from ((TJ), 6 < and c < 0. It follows that a < 1 and d < 1 
in ([9]). Furthermore, since the map T is strongly competitive, a > 0, d > 0, b < 0, c < 0. Hence, T must 
satisfy inequalities < a < 1, < d < 1 and 6c > from hypotheses ([2]) of Theorem [TJ Finally note that 
det(J T (x, y)) = ad- be > by hence 1 + (a + d) + ad- bc> 0. □ 

Denote with Q^(a,b), £ = 1,2,3,4 the four regions Qi(a,b) := {(x,y) G M 2 : a < x, 6 < y }, 
Q 2 (a,b) := {(x,y) £t 2 : x < a, 6 < y }, Qz(a,b) := {{x,y) G R 2 : x < a, y < b }, Q 4 (a,6) := {(^,y) G 
M 2 : a < x, y < b }. 

Lemma 5 Each of the sets Q\ (hi, h-z) and Qz(hi, h 2 ) contain at least one equilibrium of Eq.(LGIN). 

Proof. The existence of an equilibrium of Eq.(LGIN) in the invariant and attracting set B := [hi, h\ + 
b{\ x [h 2 , h 2 + b 2 ] from Lemma [2] is guaranteed by the Schauder Fixed Point Theorem [7] for continuous 
maps on compact, convex and invariant sets. Since B C Qi(hi, /12), such equilibrium points are elements 
of Qi(hi, h 2 ). To see that Q3 (hi, h 2 ) contains an equilibrium of Eq.(LGIN), note that the critical curves 
Ci and C2 can be given explicitly as functions of x : 

„ , s x 2 + (1 - 61 - hi) x - hi 

Ci : yi (x) = — r ,x + hi 

ci(hi - x) 



C 2 
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y2+(x) = - ( — 1 + b 2 + h 2 - c 2 x + y / (- 6 2 - h 2 + c 2 x + l) 2 + 4 (c 2 x h 2 + h 2 ) 
V2-(x) = ^ ( _1 + ^2 + hi - c 2 x - \J(- b 2 - h 2 + c 2 x + l) 2 + 4 (c 2 x /i 2 + h 2 ) 

lim yi(x) — y 2 _(x) = 00 and lim yi(x) — y 2 _(x) = — 00 . (10) 

x— *— 00 x^hi 

One can conclude from (fl~0|) and the continuity of yi(x), y 2 _(x) that there exists c < hi such that yi(c) = 
y 2 _(c). Since x = /ii is a vertical asymptote of Ci and y = h 2 is a horizontal asymptote of C 2 , it follows 
from the decreasing characters of yi(x) and y 2 _(x) that (c, yi(c)) must lie in Q3 (/ii, hi)- □ 



Then, 



Definition 1 Lei k be a positive integer. Let (xq, yo) G M 2 6e an intersection point of the graphs C\ and 
C 2 0/ two k-times differentiable functions fi, / 2 defined in a neighborhood of xq. The point (xq, yo) is a 
contact point of Ci and C 2 of order k if 



f¥\xo) = f?(x ), o<e<k 



{ f[ k \*o) + / 2 w (x ) 



(*)/ 



Note that C\ and C 2 intersect transversally at (xo, yo) if an d only if (xo, yo) is a contact point of Ci, C 2 
of order one. 
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Lemma 6 Let Po = (xq,uq) be an intersection point of two non-singular algebraic curves C\ and C2 in 
the plane that have no common component through Pq . If C\ and C2 have tangent lines at Pq that are not 
parallel to either axis, then there exists a neighborhood U of Pq for which both U D C\ and U (IC2 are the 
graphs of real analytic strictly monotonic functions y\(x) and y2(x) for x in some neighborhood of xq, and 
the order of Pq as a contact point of C\ and C2 equals the multiplicity of Pq as an intersection point of the 
algebraic curves C\ and C2. 

Proof. Without loss of generality we may assume (xo, yo) = (0, 0). Let yi(t) and y2(t) be such that UnCi 
and U fl C2 are parametrized by (i, yi(i)) and (i, 2/2^)) f° r t € I. Then y\(t) and y2{t) are real analytic on 
/. The rest of the statement follows. □ 



Lemma 7 Let p and q be real analytic strictly monotone functions defined in neighborhoods of xq and yo 
respectively, such that p(xq) = yQ, q{yo) = xq, and neither p o q nor qo p is the identity function. Let k 
be the order of (xo,yo) as a contact point of C± := {(x,y) : y = p(x)} and C2 ■= {(x,y) : x = q(y)}, and 
let d y and d x be the multiplicities of yo and xq as zeros of y — p(q(y)) and x — q(p(x)) respectively. Then 
dy — d% — k • 

Proof. Consider the function (j)(y) := y —p(p(y)), where p is the inverse of the strictly monotone function 
p. Clearly, cj)(y) = in a neighborhood of yo- It follows that (j) (yo) = for £ > 0. In particular, 4>'{yo) = 
from which we get 

_ 11 

p'{yo) = y- = Tl r ■ (11) 

p (p{yo)) p'{xo) 

Now consider the function ip(y) = y — p(q(y))- Since yo is a root of il>(y) of multiplicity d y by hypothesis, 
we must have ip^\yo) = for < I < d y and ip^ dy \yo) 7^ 0. In particular, ip'(yo) = from which we have 

q'(y ) = ^r^ = ^— v (12) 

p {q{yo)) p'{xo) 

From (fTT|) and (fT2|) . it follows that p'(yo) = q'(yo)- Similarly, one can show that p^\yo) = Q^(yo): 2 < 
£ < d y . However, since "0 (2/0) 7^ 0, we must have p^ dy \yo) 7^ Q^ dy \yo)- Hence 



P {e) (yo) = q {e) (yo), 0<£<d y 
P (dy) (yo) + q {dy) (yo) 



(13) 



It is a direct consequence of Definition Q] and (|13p that d y = k. By a similar argument using the equation 
defining C2, one can show that d x = k. □ 

The main result of this paper is the following. 

Theorem 3 The following statements are true: 

(i.) Eq.(LGIN) has at least one and at most three equilibria in [0, oo) 2 . The set of equilibrium points in 
[0, oo) 2 is linearly ordered by ^ se . 

(ii.) If Eq.(LGIN) has exactly one equilibrium in [0, oo) 2 , then it is is globally asymptotically stable. 

(iii.) //Eq.(LGIN) has three distinct equilibria in [0,oo) 2 , say (xe, y^), I = 1,...,3, with 

(xi,yi) -< se (x2,y 2 ) ^se (^3)^3), then (xi,yx) and (x3,y 3 ) are locally asymptotically stable, while 
(^2,^2) * s a saddle point. The global stable manifold 0/(^2,^2) ^ s the graph of a continuous increasing 
function of the first variable with endpoints in the boundary of [0,oo) 2 ; which is a separatrix of the 
basins of attraction of (xi,yi) and (xs,y 3 ). 
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(iv.) If there exist exactly two equilibria in [0,oo) 2 ; then one is locally asymptotically stable and the other 
is a nonhyperbolic fixed point. If (xi,yi) and (x 2 ,y 2 ) are the two equilibria with (xi,yx) ^ se (x 2 ,y 2 ), 
then Q2(x\, yi) is a subset of the basin of attraction of(xi,yi), and Q±(x 2 , y 2 ) is a subset of the basin 
of attraction of (x 2 , 2/2)- 

Proof. 

(i.) It is a consequence of Bezout's Theorem (Theorem 3.1, Chapter III in [TU]) that the hyperbolas C\ 
and C2 given in ([3|) must intersect in at most four points. Since intersection points of C± and C2 
are precisely the equilibrium points of Eq.(LGIN), it follows that Eq.(LGIN) must have at most four 
equilibrium points. Lemmas [2] and [5] guarantee that at most three of these equilibria can lie in the 
invariant attracting box B := [hi, h 2 ] X [hi + 61, h 2 + b 2 ] and hence in Q\ (hi, h 2 ), which proves the 
first part of statement (i.). To see that the set of equilibrium points in [0, oo) 2 is linearly ordered by 
^ se , note that the equilibria are precisely the intersection points of the decreasing critical curves Ci 
and C2 given in Lemma [3j 

(ii.) Let T be the map of Eq.(LGIN), and let (x,y) G [0,oo) 2 . By Lemma[2]and since T is competitive, 
(hi,h 2 + b 2 ) ^ se T(hi,h 2 + b 2 ) ^ se T(h 1 + b 1 ,h 2 ) ^ se (hi + bi, h 2 ). By induction, T n (hi, h 2 + b 2 ) < S e 
T n+1 (hi,h 2 + b 2 ) <se T n+1 (/ii + h,h 2 ) < se T n (h x + h,h 2 ) for n = 1,2, .... Then the sequences 
{T n (hi, h 2 + b 2 )} and {T n (h\+bi, h 2 )} are respectively monotonically increasing and decreasing, they 
are bounded and converge (this follows from the bounded and monotonic coordinate-wise character 
of the sequences). Since both must converge to a fixed point, the limit point is the unique fixed 
point of T. Again by LemmaO for (x, y) G [0, oo) 2 , T(x, y) G B, hence (hi, h 2 + 62) ^se T(x, y) < se 
(hi + bi,h 2 ), and by induction, T n (hi,h 2 + b 2 ) ± se T n+1 (x,y) < se T n (h x +h,h 2 ) for n = 1,2,.... 
Thus T n (x,y) — > (x,y) as n — ► 00. This proves global attractivity of (x, y) Stability follows from 
the fact that T n (h x ,h 2 + b 2 ) ^ se T n+e (x,y) < se T n (h x + bi,h 2 ) for £ = 1, 2, . . ., n = 1,2,..., and 
dist(T n (/n, h 2 + b 2 ),T n (hi + bi, h 2 )) -» 0. 

(iii.) Suppose Eq.(LGIN) has three distinct equilibria in [0, oo) 2 , say (x£, y^), I = 1, ...,3, with (xi,y^) < se 
(%2,y 2 ) ^se (x~3, Vs)- Note that the equilibria are precisely the intersection points of the critical curves 
Ci and C 2 given in ([3]). From Bezout's Theorem, Lemma[5]and Lemma[6l each equilibrium in [0, oo) 2 
must be a contact point of Ci and C 2 of order one. It follows from the remark after Definition Q] that 
Ci and C 2 must intersect transversally at each of the three equilibria in [0, oo) 2 . Furthermore, solving 
for y and x respectively in the equations defining C\ and C 2 in ([3|) gives that the vertical asymptote 
of C\ is x = hi and the horizontal asymptote of C 2 is y = h 2 . The asymptotes guarantee that in 
order to have three intersection points in [0,oo) 2 , the slopes of the functions yi(x) and y 2 (x) of C\ 
and C 2 respectively, must satisfy the relations y'i(xi) < y' 2 (xi), y[(x 2 ) > y' 2 (x 2 ) and y'x^xz) < y' 2 (^3) ■ 
Theorem [1] then gives that (xi,y t ) and (x~3,y 3 ) must be locally asymptotically stable, while (x 2 ,y 2 ) 
must be a saddle point. The rest of the proof follows from Theorem 4 and Corollary 1 in |12j . 

(iv.) Suppose Eq.(LGIN) has two distinct equilibria in [0, oo) 2 , say, (xi,y t ) and (x 2 ,y 2 ). Note that (x~i,y~i) 
and (x 2 ,y 2 ) cannot both be contact points of C\ and C 2 of order one. Indeed if they were, then as 
a result of the remark after Definition Q3 Ci and C 2 would have to intersect transversally at (x~i,y~i) 
and (x 2 ,y 2 )- Thus by Theorem [H both equilibria would have to be locally asymptotically stable. 
But this would lead to a contradiction since by Theorem 4 of [3j, at least one of the equilibria has 
to be a non-attracting equilibrium. Suppose (x~i,y~i) and (x 2 , y 2 ) are contact points of Ci and C 2 of 
orders two and one respectively. Then by the remark after Definition [H Ci and C 2 must intersect 
tangentially at (xi,yi) and transversally at (x 2 ,y 2 ). Hence by Theorem CD, Cei,I7i) must be a non- 
hyperbolic fixed point and (x 2 ,y 2 ) must be locally asymptotically stable. The statement on basins 
of attraction is a consequence of T being competitive and the fact that Q 2 (xi,yi) and Qi(x 2 ,y 2 ) 
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are invariant sets for T. Indeed, if (x,y) G Q 2 (xi,yi) then T n (x,y) € Q 2 {xi,y{). Since {T n (x,y)} 
converges to a fixed point by Theorem [21 it must converge to the only fixed point in Q 2 (xi,yi), 
namely, (xi,yi). A similar argument applies if (x,y) € Qi(x 2 ,y 2 ). 

□ 



5 A global attractivity result 

A version of the following theorem was proved first in [9]. It is given here for easy reference. 

Theorem 4 Let [a, b] and [c, d] be intervals of real numbers and let f : [a, b] — > [a, b] and g : [c, d] — > [c, d] 
be continuous functions that satisfy the following properties (a) and (b) ; 

(a) f(x,y) is non- decreasing in x and non-increasing in y, and g(x,y) is non-increasing in y and non- 
decreasing in x. 



(b) For every (m,M), (m,M) € [a, b] x [c,d], 



m = f(m,M), m = f(m,M), m = m 

M = f(m,M), M = g(m,M) M = M 



(14) 



Then the system of difference equations 

I {t 1 ^] n = 0,1,2,..., (x , yo )e[a,b)x{c,d) (15) 

yn+l — y\Xmyn) 

has a unique equilibrium (x,y) in [a,b] x [c,d], and the unique equilibrium is a global attractor. 
Next we give a sufficient condition for the existence of a unique equilibrium. 
Theorem 5 // at least one of the following conditions is satisfied 

(a) 1 — b\ + hi + ci /i2 > and 1 — 62 + ^2 + C2 h\ > 

(16) 

(b) cic 2 <l 

then Eq.(LGIN) has a unique equilibrium in [0, oo) 2 . 
Proof, 
(a) Suppose 

1 - bx + hi + ci h 2 > and 1 - b 2 + h 2 + c 2 hi > 0. (17) 

From direct inspection of ([7]) one can see that the functions f(x, y) and g(x, y) in T(x, y) = (f(x, y),g(x, y)) 
satisfy hypothesis (a) of Theorem HI To verify hypothesis (b) of Theorem [H note that for T as in 
Q, the system of equations in (fT4"|) is given by 



(18) 



bim , ^ _ b 2 m 

m = = + /ii m = 7-—— — + h 2 

1 + m + ciM l + m + c 2 M 

> and 

biM ; — b 2 M 

M = -4 — + hi \ M = = h h 2 

1 + M + cim ) 1 + M + c 2 m 
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Algebraic manipulation of the equations in (|18p yields the equation 

c 2 (M_- m) [(m - /ii) + (M_- fa) + (1 - 61 + *n + a h 2 )] 

+ ci(M - m) [(fn - h 2 ) + (M - h 2 ) + (1 - b 2 + h 2 + c 2 hi)\ =0 [ ' 

Using hypothesis (|17|) and the facts m, M > h\ and fn, M > h 2 , we obtain that in equation (|19|) . 
m = M and m = M which shows that hypothesis (b) holds. Then Theorem H] implies that there is a 
unique equilibrium. 

(b) Suppose ci c 2 < 1. Applying the transformation x = X + hi, y = Y + h 2 to the equations defining 
the sets Ci , C 2 in Lemma [3] gives rise to the equations 



(20) 



(a) (l- Cl c 2 )X 4 + BX 3 + CX 2 + DX + b\h\ = 

(b) (l- Cl c 2 )Y 4 + BY 3 + CY 2 + + b\h\ =0 

where -B, C, D, B, C and D are functions of 61, b 2 , ci, c 2 , di, d 2 . We consider two cases. 

Case 1 : Suppose c\c 2 <1. Since the leading coefficient of ([20a ) is 1— c\ c 2 and the constant coefficient 

b 2 h 2 

is b\ h 2 , the product of the roots of (f20h ) is - — - — - — > 0. Similarly, the product of the roots of ([20b ) 

1 - C\ c 2 

b 2 h 2 

is — - — - — > 0. From Lemma [5l each equation in (I20p must have a positive root a and a negative 
1 - ci c 2 

root /3. Furthermore, none of the two equations can have complex roots. Indeed if, say, (|20r ) had a 
pair of complex conjugate roots z and z, then the product of the roots of ([20h ) given by a/3|z| 2 would 
be negative, which is impossible. It follows that the set of roots of each equation in (I20j) must have 
the sign structure {+, +, — , — }. As a result, the set of equilibria of Eq.(LGIN) must have one of 
the sign structures (i.) {(+, +), (-, -), (+, +), (-, -)} or (ii.) {(+, +), (-, -), (+, -), (-, +)}. 
The sign structure given by (i.) is not possible since otherwise, there must exist two equilibria of 
Eq.(LGIN) in and two equilibria outside . By the proof of Theorem [3] part iv., one of the 
equilibria in must have multiplicity two while the other equilibrium in must have multiplicity 
one. Thus the sum total of the multiplicities of all the equilibrium points must be at least five, 
contradicting Bezout's Theorem. Hence when c\ c 2 < 1, the only possible sign structure for the set 
of equilibria of Eq.(LGIN) is (ii.), which implies that Eq. (|20p has a unique equilibrium in Qi(0, 0), 
and hence Eq.(LGIN) has a unique equilibrium in Qi{h\,h 2 ). 

Case 2 : Suppose c\c 2 = 1 : In this case, the equations in (i20"j) reduce to 

{l-b 1 -c 1 + b 2 c 1 )X 3 + C1X 2 + D 1 X + b\h\ = 

_ _ ( 21 ) 
- ci (1 - h - ci + b 2 a) Y 3 + Ci Y 2 + Di Y + b\ h\c\ = 

where C\, D\, C\ and D\ are functions of by, ci, hi, b%, c 2 , h 2 . We consider three cases: 1 — t»i — 
c\ + b 2 c\ is negative, zero, and positive. If 1 — b\ — c\ + b 2 c\ 7^ 0, the product of the roots of the 

b 2 h 2 

first equation in (I2ip is ; — and the product of the roots of the second equation 

1 — bi — c\ + o 2 ci 

b 2 h 2 c\ 

in (1211) is ; — - — : . Note that from Lemma [51 we have that each equation in (12111 must 

1 - 61 - ci + b 2 a H """" 

have a positive root and a negative root. Thus the equations in (|2ip cannot have complex roots 
due to their cubic nature, since complex roots of real polynomials always occur in conjugate pairs. 
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If 1 — bi — c\ + 62 c\ < 0, the sets of roots of the equations in ([2T|) must have sign structures 
{+, — , — } and {+, +, — }, respectively, by an argument similar to the one already used in case 1. 
Hence the set of equilibria of Eq.(LGIN) must have the sign structure {(+, +), (— , — ), (— , +)}. If 
1 — b\ — c\ + 62 ci > 0, the sets of roots of the equations in (f2Tj) must have sign structures {+, +, — } 
and {+, — , — }, respectively. Hence the set of equilibria of Eq.(LGIN) must have the sign structure 
{(+, +), (— , — ), (+, — )}. If 1 — 61 — ci + 62 ci = 0. the equations in (|2"T|) reduce to a pair of quadratic 
equations. It follows from Lemma [5] that each equation must have a positive root and a negative 
root. So the set of roots of each equation must have the sign structure {+, — }. Hence the set of 
equilibria of Eq.(LGIN) must have the sign structure {(+, +), (— , — )}, which implies that Eq. ([2"Uj) 
has a unique equilibrium in (Ji(0, 0), and hence Eq.(LGIN) has a unique equilibrium in Q\{h\,h2). 

□ 
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